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■ ABSTRACT 
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, We investigate the effect of line of sight temperature variations and noise on 

■ two commonly used methods to determine dust properties from dust continuum 
O: observations of dense cores. One method employs a direct fit to a modified black- 
i-C ■ body SED; the other involves a comparison of flux ratios to an analytical pre- 
diction. Fitting fluxes near the SED peak produces inaccurate temperature and 
dust spectral index estimates due to the line of sight temperature (and density) 
variations. Longer wavelength fluxes in the Rayleigh- Jeans part of the spectrum 

'— ( £ 600 /iin for typical cores) may more accurately recover the spectral index, but 
y—i . both methods are very sensitive to noise. The temperature estimate approaches 

^ . the density weighted temperature, or "column temperature," of the source as 

. short wavelength fluxes are excluded. An inverse temperature - spectral index 

^ . correlation naturally results from SED fitting, due to the inaccurate isothermal 

(sj . assumption, as well as noise uncertainties. We show that above some "threshold" 

temperature, the temperatures estimated through the flux ratio method can be 
O ■ highly inaccurate. In general, observations with widely separated wavelengths, 

£> . and including shorter wavelengths, result in higher threshold temperatures; such 

•i-H . 

^ , observations thus allow for more accurate temperature estimates of sources with 

temperatures less than the threshold temperature. When only three fluxes are 
available, a constrained fit, where the spectral index is fixed, produces less scatter 
in the temperature estimate when compared to the estimate from the flux ratio 
method. 
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Introduction 



Some of the coolest regions in molecular clouds are dense, starless cores. These dust en- 
shrou ded objects are often in the process of forming one (or a few) protostar(s) (IBenson fc Myers 
19891 ). Determining the physical properties of cores, such as temperature, composition, and 
density, is necessary for a complete understanding of the environmental conditions prior to 
the formation of a star or protostar. There has been much progress in the study of cores 
containing central protostars, including the success of theory in explaining t he variety of 



emergent spectral energy distributi ons (SED) as an evolutionary sequence (e.g. lAdams et al. 



19871 : lLadalll987l : lAndre et al.lll993l ). On the other hand, the structure and evolution of cores 
that have yet to form a central protostar is not as well understood, and is thus an active 
area in current star formation research. 

Dust presents one avenue to observationally investigate starless cores. Dust is prevalent 
in the ISM, and it is responsible for most of the extinction of starlight. For cores positioned 



, „„„ ^ 

(Lada et al. 


1994; 


Alves et al. 


2001) 



scattered light from co res surrounded by diffuse background radiation may be used to deter- 



mine the dust content (IFoster fc Goodmanll2006l ). Dust can also be directly detected through 
its thermal emission. Since the temperatures of the cores are ^ 15 K, the emergent contin- 
uum SED peaks in the far- infrared (FIR) or sub-millimeter wavelength regimes (see Fig. [TJ. 
Ground and space based observations by SCUBA, MAMBO, Bolocam, 2MASS, IRAS, ISO 
and Spitzer have detected dust emission fr om many environments, and they have provided 



much information abou t starless cores (e.g. I Ward- Thompson et al.ll2002l ; ISchnee et al.ll2007 



Kauffmann et al.ll2008l ). The upcoming Planck and Herschel missions, which are capable of 
FIR observations, are also well suited for detecting dust emission. Thus, a thorough con- 
sideration of the nature of dust continuum emission, and the uncertainty associated with 
measuring it, is timely. 

The main characteristics of the dust that determine the form of the emergent SED 
are the column density, temperature, and emissivity. Observationally quantifying these 
characteristics should constrain models of dense starless cores. Radial densi t y profiles ar e 
often compared with a stable isothermal Bonnor-Ebert sphere ( IBonnorlll956l : lEbertlll955l ). 
for which the volume density is constant near the c enter, but drops as r~ 2 at larger radii 
(e.g. iBacmann et al.ll2000l ; ISchnee &: Goodman!l2005l ). The density gradients may vary from 
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core to core, which in turn may (or may not) be due to an evolutionary sequence as cores 
continually collapse to form a protostar. 

Though gas temperatures in a stable Bonnor-Ebert sphere are constant, theoretical 
results have s uggested that dust temperatures decrease towards the center to values as low 
as ~ 7 K (e.g. lLeungill975l ; lEvans et al.ll200ll ; IZucconi et al.ll200ll ). And, recent observational 
inves tigations have in f a ct identified cores with such gradients in the dust temperatures 



(e.g. 



Schnee et al.l 120071 ; IWard-Thompson et al.l [2002]) . Though the dust mass is only a 
fraction of the gas mass (~ 1/100), gas temperatures may also exhibit gradients, due to th e 



coupling between dust and gas at high enough densities (lGoldsmithll200ll ; ICrapsi et al.ll2007l ). 
The (recent and upcoming) availability of higher quality observational data will require a 
thorough interpretation of emergent SEDs to accurately assess the temperature, as well as 
density, profiles of the observed sources. 

The common assumption is that the emergent SED from interstellar dust is similar to 
the Planck functio n of a blackbody, modified by a power-law dependence on the frequency 
(IHildebrandl Il983l ) . The spectral index of the dust emissivity power-l aw, /3, is dependen t 
on the bulk and surface properties of the dust grains. As shown by iKeene et al.l (119801 ) . 
observations limited by sparse flux sampling may be consistent with various SEDs described 
by different values of (3. A precise estimate of the value of (3 is necessary to accurately derive 
other properties of the observed source, such as the temperature and the mass of a cold core. 
The emissivity-modifie d blackbody spectrum is the basis for m a ny analyses of dust properties 



in observed core s (e.g Kramer et al. 



Kirk et al.l 120071 ). 



20031 : ISchnee et al.ll2005l : IWard-Thompson et all 12002 



Dupac and coworkers fit observed FIR and sub-millimeter fluxes with a modified black- 
body spectrum, and they suggested that j3 decreases with inc reasing temperatures, fr om ~ 2 



in cold regions to 0.8 - 1.6 in warmer regions (T ~ 35 - 80 K; lDupac et al.ll200ll . 120021 . 120031 ). 
Such analyses may be sensitive to the simplified assumption of a constant te mperature along 
the lin e of sight. Using radiative transfer calculations of embedded sources, iDoty fc Leung 
(119941 ) demonstrated that the accuracy of the parameters estimated from measured fluxes 
is sensitive to the precise nature of the source (e.g. opacity, temperature distribution); they 
found that the Rayleigh- Jeans (R- J) regime of the emergent spectrum is better suited for an 
accurate determination of the dust spectral index. IDoty fc Palottil (120021 ) found that the use 
of flux ratios to estimate the spectral index is sensitive to which wavelengths (of the given 
fluxes) are used in the ratio; they also found that is more acc urately determined when 
fluxes at longer wavelengths are used in a fit. I Schnee et al.l (120061 ) also showed that various 
ratios of fluxes (with different wavelengths) give different estimates for dust temperature and 
column density, due to an inaccurate isothermal assumption. 
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Here, we systematically investigate how line of sight density and temperature varia- 
tions, similar to those in dense cores, as well as noise uncertainties, affect the temperature 
and spectral index estimated from IR and sub-millimeter observations. We focus on two 
commonly employed methods. The first method uses a direct fit of a modified blackbody 
SED. For the second method, ratios of the fluxes are used to determine the temperature and 
(3 from an analytical prediction. Both methods usually rely upon an assumption of constant 
temperature along the line of sight. Using simple radiative transfer calculations of model 
sources, and Monte Carlo experiments, we assess how well the resulting temperature and 
spectral index estimates recover properties of known input sources. 

This paper is organized as follows. In the next section (§2), we present the analytical 
expression of the power-law modified blackbody spectrum, and briefly introduce two well 
known methods used to estimate the dust properties from IR and sub-millimeter continuum 
observations. In our analysis of the two methods, we consider numerous scenarios typical 
of observations of star forming regions. Table [1] shows the particular scenario considered 
in each subsection, and may be used as a brief guide to §3 - §6. We begin our analysis 
by considering non-isothermal sources in the ideal limit where a large range of fluxes at 
different wavelengths are available for fitting an SED. We then systematically exclude fluxes, 
culminating with the scenario where only a few fluxes are available, in which case the flux 
ratio method is employed. Throughout our analysis, we also consider the effect of noise 
in the observations, of both isothermal and non-isothermal sources. In §[3j we describe the 
method to estimate source temperatures using direct SED fitting; we investigate how line 
of sight variations, noise, and the sampling of different regions of the emergent SED affect 
the resulting temperature estimates. We also discuss our findings in the context of recent 
published works. In the following section (§3]) we analyze the flux ratio method focusing 
on the effect of noise, through the use Monte Carlo simulations. We then compare the two 
methods using a radiative transfer simulation to model the emission from an isolated starless 
core in §5J After a discussion in §H1 we summarize our findings in §7J 

2. Dust Emission: Common Assumptions and Methods 

2.1. Isothermal Sources 

The emergent continuum SED due to dust is often expressed analytically as the product 
of a blackbody spectrum B U (T) at the dust temperature T and the frequency dependent dust 
opacity k v . The observed flux density associated with this SED takes the form 



S v = QB u (T)k u N, 



(1) 
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where f2 is the solid angle of the observing beam, and N is the column density of the emitting 
material. The opacity k„ is empirically determined to have a power law dependence on the 
frequency (jHildebrandlll983l ): 



Kb 



v 

K>0 | — 

.^0 



(2) 



The spectral index (3 depends on the physical and chemical properties of the dust. For 
silica te and graphite dust composition common in much of the ISM, (3 ~ 2 (IDraine fc Lee 
19841 ) . However, observations have sho wn that j3 can reach values as low as ^ 1 and as hig h 
as ^ 3 in various environments (e.g. lOldham et al.lll994c iKuan et al.lll996l ; lMathislll990l ). 
Indeed, the spectral index is a key parameter and accurately determining its value, along 
with the column density iV and the temperature T, is crucial for a thorough description of 
dust properties in an observed region. Those are the three parameters that are required to 
accurately describe an observed flux density (per beam, i.e. S u /Q). 

Figure [TJ shows SEDs from a 20 K source with different values of (3, but constant column 
densities. For comparison, the SED from a 5 K source with the same column density and (3 = 
2 is also shown. The SEDs are all calibrated using an equivalent k = ^23ogh Z) which is why 
all the 20 K SEDs intersect at 230 GHz. Sources with higher spectral indices (f3) produce 
SEDs with steeper slopes at long wavelengths (in the R-J regime), increasing peak fluxes, 
and shorter peak wavelengths. Similar to Wien's Law for a pure blackbody, for a given value 
of (3, a modified Wien's Law indicating the wavelength corresponding to the peak in S u , 
Amax, can be determine d numericallylj] For (3 = 2, A max ~ (2900/zmi\~)/T, and for (3—1, 



X max ~ (3670 fan #)/T. botv fc Palott] d2002h find that X max = (4620e-°- 2357/3 fimK)/T is 
a good fit for 1 < (3 < 2. 



2.2. Non-Isothermal Sources: A Simple Example 

Equations ([I])-© describe the spectrum emitted from dust at a single temperature T. 
For a 3D source with various dust characteristics, the emergent SED will be a combination 
of SEDs from all the dust in the source. For optically thin emission, the emergent SED is 
simply the integrated SED from each dust grain. Here, we briefly consider the effect of using 
equations ([I])-® to characterize an SED from a source with two different dust populations. 



1 This modified Wien's Law gives the wavelength that corresponds to the peak in S v . In other texts, A. 
sometimes refers to the wavelength corresponding to the peak of S\. 
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This relatively simple analysis is a prelude to the effect of line of sight temperature variations 
on commonly employed methods to estimate dust properties. 

Figure [2] shows the emergent SED from a source with two populations of dust grains, 
along with the SED from each individual component: the temperature and column density 
of the cool component is T\ = 10 K and Nx, respectively; for the warm component, T2 = 15 
K and N2 = O.lAq. Physically, such a system is similar to a 10 K dense core surrounded by 
a 15 diffuse envelope. The spectral indices for both components are set to (3 = 2. 

The peak of the emergent SED in Figure [2] occurs at X max = 251 /xm. Using the 
modified Wien's Law for (3=2, the temperature of an isothermal source that would produce 
a spectrum which peaks at that wavelength is 11.6 K. Though this temperature occurs 
somewhere between the temperatures of the two isothermal sources that contribute to the 
emergent SED, in practice knowledge of the peak of the SED of an unknown source is not 
easily determined. Further, it is not obvious how one should interpret the temperature 
assigned to a source which itself is not isothermal. 

In subsequent sections, we assess how well emergent SEDs are described by equations 
CO)-©, and what information about the source temperature can be garnered from continuum 
observations that span different regions of the SED. We also consider the more realistic 
constraint of limited sampling of the emergent SED. One question we aim to address, for 
instance, is whether fluxes in different parts of an emergent SED, such as the Wien or R-J 
regimes, are preferable for determining the dust properties. 

Two commonly employed methods to determine the properties of dust from continuum 
observations are: (1) a direct fitting of equations fll} & ©; and (2) the use of ratios of 
observed flux densities at 2 or more wavelengths. For isothermal sources, and with ideal 
observations with no uncertainties, both methods will accurately recover T, (3, N, and k v . 
However, all sources are unlikely to be isothermal, and even the most accurate observations 
include some level of intrinsic noise. In the following sections, we quantify how these factors 
affect the accuracy of the derived parameters. After describing the two methods in further 
detail, we use simple numerical experiments to evaluate the accuracy of the methods in 
determining the dust properties from observations of star forming cores. 

3. Direct SED Fitting 

A minimized \ 2 fit of equations ([I])-© to a number of observed fluxes can be performed 
to estimate the dust properties. There are essentially three parameters to be fit, the tem- 
perature T, the spectral index (3 and the absolute scaling, which is just the product of the 
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column density N and the opacity at a given frequency k,q. Since a fit will only produce the 
scaling (which is the optical depth at a particular frequency, e.g. r 2 3o = -/Vk 2 3o at 230 GHz), 
other assumptions and/or techniques are necessary to obtain estimates of N and k,q. For 
example, if the opacity at a wavelength is known (e.g. k = k 2 3o for v = 230 GHz), then 
the fit can estimate N directly. Extinction studies are another avenue to estimate N; the 
level of attenuation (usually from optical and NIR observations) due to dust in dark clouds 
in front of the st ellar background is directly related to the column density of the dust (e.g. 



Lada et al.lll994j ). This method is advantageous since it provides an independent estimate of 
N, but also requ i res assumptio ns, such as the ratio of total-to-selective extinction Ry (e.g. 



Hildebrandlll983l ; lMathislll990l . and references therein). Further, to obtain the total column 



density along the line of sight, and not just that of the dust, an additional assumption of the 
dust-to-gas ratio is required. In our analysis, we will assume that only IR and sub-millimeter 
observations are available, and thus will limit our analysis to the estimation of T and /3. Our 
focus here is to investigate how well a given method can reproduce temperatures and spectral 
indices only; estimation of the absolute column density and opacity is beyond the scope of 
this work. 



3.1. Effect of Line of Sight Temperature Variations 

3.1.1. Two-Component Sources 

We begin by considering ideal (i.e. error-free) observations of simple two-component 
so urces. Fitting experi ments involving sources with two dust populations have been explored 



by iDupac et al.l (120021 ); their aim was to determine the amount of cold dust, along lines of 
sight with warmer dust, that is necessary to reproduce the fit results of their observations. 
Here, we simply evaluate the resulting fits when fluxes at different wavelengths are available. 

The temperatures of the cold and warm media are T\ and T 2 , respectively, and the 
column density ratio is NijN\. Such systems are analogous to isothermal (spherical) dense 
cores surrounded by warmer envelopes. Modified blackbody SEDs (eqns. [I]- [2]) are con- 
structed for the two media, both with /3 — 2. A particular example SED of this general case 
is shown in Figure [2J We then fit equation (pQ) to the integrated SEDs, solving for T and (3 
(as well as the scaling factor Nkq). Figure [3] shows the fit temperatures and spectral indices 
from observations of a variety of two-component systems, along with the wavelength range 
of the fluxes considered in the fit. 



Since T varies along the line of sight, the best fit T will likely not be equal to the tem- 
perature of one of the two sources. One characteristic temperature of this two-component 



- 8- 



medium is the density weighted temperature (e.g. iDoty fc Palottil 120021 ) . Since this density 
weighted temperature is analogous to the column density, we will call it the "column temper- 
ature," T co i. The estimated temperature from a fit can be compared with this true column 
temperature. 

As shown in Figure [3j the best fit temperature is systematically too high when all fluxes 
at (integer) wavelengths between 10 - 3000 fim are considered in the fit. In fact, when 
the temperature difference between the two components is large, the best fit temperature is 
actually larger than the warmer medium, as in "2COMPc" and "2COMPd," lines of sight 
containing dust at 10 and 20 K. Further, when all wavelengths are considered, the fit value 
of /3 is always lower than the actual value of 2. However, when shorter wavelengths are 
systematically excluded from the fit, the best fit temperature decreases and approaches the 
column temperature. The best fit f3 also approaches the model value of 2 when fluxes in 
the R-J wavelength r egime are the only ones used in the fit, consistent with the findings of 
Doty fc Leung (11994I ). In practice, for deriving dust properties from the emergent SED, it 
may be necessary to exclude fluxes with A ^ 100 Atm due to the contrib ution of embedded 
sources as well as transiently heated very small grains (]Li &: Draindl200ll ). depending on the 
environment. 



3.1.2. Cores with Density and Temperature Gradients 



Recent theoretical and observational studies have indicated that the d ust temperature in 
starless cores decreases toward the center, reaching low values £j 7 K (e.g. lEvans et al.ll2001 



Crapsi et al.ll2007l ; ISchnee et al.ll2007l ; I Ward- Thompson et al.ll2002l ). We thus investigate the 
emergent SED from cores containing temperature gradients like those observed, and whether 
any useful information can be obtained from fitting a single power-law modified blackbody 
spectrum to that SED. 

To constru c t the m odel cores, we use the density and dust temperature profiles presented 
by lEvans et al.l (120011 ). who performed radiative transfer simulations on a variety of model 
cores with a range of density profiles. Though the resulting dust temperature profiles T(r) 
are sensitive to the model density profiles, the relationship between T and column density 
N is relatively unifo rm between models with different (volume) density profiles (see Fig. 9 
of lEvans et al.ll200ll ). In our analysis of emergent SEDs from starless cores, we construct 
two cores, with temperatures ranging between 8 - 12 K and 5 - 12 K; the column densities 
are as indicated in Figure 9 of lEvans et al.l (120011 ): N increases with decreasing temperature 
(and thus with decreasing core radius). At the outer edge of the core, at a temperature of 
12 K, the column density is set to 2xl0 21 cm -2 ; the temperature drops to 8 K in model 
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Core 1 and to 5 K in model Core 2, with column densities of 1.25 xlO 22 cm~ 2 and lxlO 23 
cm -2 , respectively. Besides these isolated cores, without any surrounding medium, we also 
consider cases where the cores are surrounded by an envelope with a temperature of 20 K 
and a column density of lxlO 21 cm -2 . 

As described in §3.1.11 we begin by fitting the emergent SED assuming that fluxes at 
various wavelength ranges are available. Figure 0] shows the resulting best fit temperatures 
and spectral indices for the two cores. The SEDs from the cores without an envelope are 
analogous to an SED obtained by accurately subtracting off flux due to larger scale emission 
from the surrounding region, or an SED from a truly isolated core. When equation ([1]) is fit 
to fluxes at 100 - 600 fim, the best fit temperature is ~ 3 K (~ 15%) off from the column 
temperature of Core 1, but differs from the column temperature of Core 2 by ~ 6 K (~ 
50%). The best fit spectral index also shows large variation between the two cores. For 
Core 1 (T G 8 - 12 K) the fit (3 of 1.65 is within 20% of the model value of 2. However, 
for Core 2 (T G 5 — 12 K), the best fit (3 of 0.81 is erroneous by over a factor of 2. As 
more short wavelength fluxes (in the Wien regime) are excluded, the fits recover the model 
spectral index more accurately; the temperature estimate also decreases, approaching the 
column temperature of the core. 

When the core is surrounded by a warmer envelope, or when the flux from extended 
regions has not been properly accounted for, then the discrepancy between the best fit 
parameters and the core properties increases, as expected. For Core 2, a fit to the fluxes 
at 100 - 600 /im results in an estimate for (3 with an unphysical sign (-0.3). Including the 
envelope, the discrepancy between the fit T and the column temperature at short wavelength 
fluxes increases by ~50%, compared with the cores without the envelope. 

We performed such fits for cores surrounded by more diffuse envelopes, with a column 
density that is a factor ~ 10 - 100 times lower than that of the core (lxlO 20 cm -2 ). Such 
envelopes only have a slight effect on the fit temperatures, because the column temperatures 
are not significantly different compared to the isolated cores. The envelope does have an 
appreciable effect on the best fit (3 when fluxes near the peak of the SED are considered 
in the fit. However, excluding short wavelength fluxes still recovers the true spectral index 
reasonably accurately, as Figure [T] would suggest. 

Figure [5] shows the emergent SED from Core 2 (T G 5-12 K) without an envelope, 
along with the results from two fits, one to the fluxes from 100 - 600 /im, and the other to 
fluxes from 1-3 mm. The fit at shorter wavelengths does indeed reproduce the peak and 
shorter wavelength fluxes of the emergent SED reasonably well, but severely overestimates 
the fluxes at longer wavelengths. On the other hand, the fit to long wavelength observation 
reproduces the R-J tail of the spectrum accurately, but underestimates the fluxes at all 
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shorter wavelengths. 

The peak of the emergent SED from Core 2 occurs at A = 324 /xm. From the mod- 
ified Wien's Law (for (3=2) \ max = (2900 fj,mK)/T, the temperature associated with this 
wavelength is 8.9 K. This value is closer to the maximum temperature, 12 K, of the source, 
than the column temperature T co ;=6.2 K. In this case, T co i is dominated by the high density, 
cold center, whereas the peak of the SED is primarily influenced by the warmer, low density 
regions of the core. Since the SED has an exponential dependence on the temperature in 
the Wien regime, which, for such a cold core, goes up to ~ 100 jttm, even low density regions 
may dominate the total SED emerging from a region, due to the higher temperatures. 



3.1.3. Summary of the Effect of Line of Sight Temperature Variations 

As expected, the emergent SED at wavelengths near the peak is poorly described by a 
simple power-law- modified-blackbody, due to te mperature variations a l ong the line of sight , 



as previously documented in the literature (e.g. iDoty fc Palottil 120021 ; ISchnee et all 120061 ). 
Nevertheless, the fits can still reveal properties of the observed source, namely that the 
best fit temperature approaches the column temperature and the best fit f3 approaches the 
model value as shorter wavelengths are excluded. There are also other revealing trends that 
warrant further investigation. First, the systematic exclusion of short wavelength fluxes 
results in different T and (3 estimates; for an isothermal source, the fit T would always 
be the same (and equal to the temperature of the source). If this trend also occurs when 
there are only a few observations, then short wavelength observations can still be used to 
determine whether a source is isothermal or not. Second, there also appears to be an inverse 
T - (3 relationship: whenever T is overestimated, (3 is underestimated. A similar relationship 



(Dupac et al. 


2001. 


2002, 


2003) 



al. would be warranted to rule out that the inferred anti-correlation is simply due to line of 
sight effects. As we discuss in the next section ( §3.2p . an inverse T - (3 trend may also arise 



from SED fits due solely to noise in the observations. 



3.2. Effect of Noise 

Uncertainties in observed fluxes may also lead to incorrect temperature and spectral 
index estimates from SED fitting. To assess the effect of noise, we first consider fluxes from 
isothermal sources with modest 5% uncertainties in each observed flux. A number fluxes 
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from a range of (integer) wavelengths are considered in the fit: 100 - 600 /mi, 500 - 1000 /mi, 
and 1000 - 1500 /mi. In these Monte-Carlo experiments, each flux is modified by a random 
value drawn from a Gaussian distribution, with a = 0.05. 

Figure [6] shows the best fit T and (3 estimates for isothermal sources, one with T = 10 
K and (3=2, and the other with T = 20 K and (3=2. Each set of noisy fluxes, spanning the 
different wavelength regimes, is generated 100 times; a modified blackbody is then fit to each 
set. As expected, there is a spread in the best fit T and (3. 

Fits from both the 10 K and 20 K source show little scatter when only fluxes with 
wavelengths 100 - 600 /im are considered in the fit, suggesting that an SED is not very 
sensitive to noise in the Wien regime. However, at longer wavelengths, there is a large 
spread in the estimated T and (3. Including fluxes between 500 - 1000 /im, the fits from the 
10 K source give (3 G 1.6 - 2.3 (within ~20% of the source value) and T G 8 - 14 (within 
~40%). This range increases when including fluxes between 1000 - 1500 /im: (3 G 1.5 — 2.6 
(within ~30%) and T G 6 — 30 (within only ~200%). That the longer wavelength fits show 
more scatter is not unexpected given the shape of an SED (see Fig. [I]). 

In the R-J tail, the SEDs from sources with different temperatures are similar in shape, 
since the slope is determined primarily by (3 (see Fig. [1]). Thus, small errors in the observed 
fluxes may result in inaccurate (3 (and thus T) fits. At shorter wavelengths, T largely 
determines the shape of the SED, so small uncertainties may be insufficient to significantly 
alter the temperature that best matches the observation in a fit. For the 20 K source in 
Figure El there is less of a difference between the spread in (3 and T for the fits to 500 - 1000 
/im and 1000 - 1500 /im fluxes. This occurs because at 20 K, the range spanning the shorter 
wavelengths (500 - 1000 /mi) is already well enough into the R-J region of the spectrum (see 
Fig. [I]), so the fit is already very sensitive to noise. 

The clear inverse (3 - T trend that emerges from SED fits to noisy fluxes is similar to 
the trends found from fits to noise-free fluxes from non-isothermal sources (as suggested by 
Figs. [3]and|4]). Apparently, whenever a fit underestimates the temperatures, the spectral 
index is overestimated, and vice versa. This effect is amplified when noisy fluxes in the R-J 
regime of the spectrum are considered in the fit, and can be understood when comparing 
SEDs with different values of (3. Consider a 20 K isothermal source, with (3 = 2, shown in 
Figure [TJ Assume the source is observed at various wavelengths, primarily at A > 1 mm, 
but that the peak is also sampled. Further, assume the noise level in the fluxes is such that 
the least squares fit preferentially obtains a (3 of 1. The peak of the SED with (3 = 1 in 
Figured] occurs at longer wavelength than the peak of the (3 = 2 SED, and so a fit T = 20 K, 
with (3 = 1, would not reproduce the peak of the observed SED well. In order for the fit to 
reproduce the peak of the SED from the 20 K source, but with a best fit /3 of 1, the best fit 
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T must be held at a larger value than 20 K. In general, therefore, when a fit underestimates 
(3, T is overestimated. 

We have shown that an inverse correlation between T and (3 can occur due to an incorrect 
assumption of isothermality, or due to intrinsic noise in the observations. Such a trend 
would of course also occur for noisy observations of non-isothermal sources. We discuss the 
combination of noise and line of sight temperature variations when we discuss the additional 
limitation of including only a small number of fluxes in §3.41 

3.3. Estimating T and /3 with Sparse Wavelength Coverages 

The fitting experiments indicate that, for lines of sight with starless-core-like tempera- 
ture and density gradients, the resulting fits produce lower temperatures and higher spectral 
indices as shorter wavelength fluxes are excluded in the fit. In all our tests so far, we have 
included all wavelengths within a given range. In practice, however, obtaining only a (small) 
number of observations at different wavelengths of a source is typically feasible. Current 
(e.g.Spitzer, SCUBA, Bolocam, MAMBO) and near future observations (e.g. Herschel and 
Planck) can provide fluxes at a number of FIR and sub-millimeter wavelengths. However, 
many of those wavebands provide fluxes at or near the peak of the integrated SEDs of dense 
cores. Thus, SED fitting may not provide accurate estimates of the (column) temperature 
and spectral index. Yet, determining whether or not a core contains temperature variations 
is instructive in itself. Further, the temperature obtained from a fit to long wavelength fluxes 
can be deemed the upper limit of the coldest region within the core. We thus investigate 
whether the identification of temperature variations along the line of sight, as well as an 
accurate estimate of the temperature limit, are still feasible when a small number of fluxes, 
including those with wavelengths near the peak of the SED, is used in the fit. We begin by 
describing the effect of fitting an SED to a limited number of noise- free fluxes. After that, 
we consider the effect of noise in observations of sources with temperature variations, in §3.41 

We construct another optically thin 2 component medium, with T\ = 10 K, T 2 = 15 
K, N 2 /N 1 = 0.02, and (3=2, similar to "2COMPb" in Figure El but with a greater density 
contrast. We perform a fit assuming fluxes were obtained at 60, 100, 200, 260, 360, and 580 
/im. The reason we choose these particular parameters, which are realistic for innermost and 
outermost regions of a core, is to assess whether any patterns emerge from fitting different 
fluxes from a source with only slight temperature variations. We choose fluxes at wavelengths 
near the peak for more direct comparisons with real observations of cores. The column of 
panels in Figure [7(a) shows the fitting results. The bottom panel of Figure [2(a) indicates the 
observed wavelengths, along with the emergent SED from the 2 component source. There are 
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no significant differences in the fitting results when excluding the 60 /im flux. However, when 
the additional flux at 100 /im is excluded, the best fit temperature decreases and spectral 
index increases. We were unable to obtain a good fit with only three fluxes; we have found 
that a minimum of four wavelengths is necessary to obtain a fit with three free parameters! 

Excluding short wavelength fluxes, even with a limited number of observations in the 
Wien regime, results in a lower estimate of T and a higher estimate of (3. We found the same 
trend regardless of what wavelengths were considered in the fit, and for a variety of systems 
with an inverse relationship between the column density and temperature. 



3.4. Implications to Recent Observations 

We now perform a similar tests to publ ished observat i ons o f cold star forming regions 



beginning with the observations presented by lStepnik et al.l (120031 ) . They analyzed a filament 



in the Taurus molecular cloud using 60, 100, 200, 260, 360, and 580 /im fluxes observed from 
IRAS and the balloon borne experiment PRONAOS/SPM. After subtracting off emission 
attributed to the surrounding envelope, they fit the 6 data points to obtain an estimate of 
the temperature and spectral index of the filament. We carry out the exact same procedure, 
but also perform fits excluding the short wavelength fluxes. The fitting results are shown in 
Figure [7(b). As with the two-component model shown in Figure [7J(a), the exclusion of the 60 
fim flux does not alter the fit. However, additionally excluding the 100 //m flux decreases the 
fit T by ~ 1 K, and increases (3 from 1.98 to 2.13. This variation is rather similar to the simple 
2 component fit, indicating that there might be a temperature variation within the filament, 
and that the actual value of the spectral index is greater than 2.1 (if the spectral index itself 
is constant in the filament). A value of 11.13 ± 1.29 K can be assigned as an upper limit 
for the column temperature of the filament. The interpretation of a temperature variation 
within the filament cannot be definitive, however, due to the large uncertainties. Within 
the uncertainties, the fit T and (3 are constant regardless of which fluxes are considered in 
the fit. More observations, preferably at longer wavelengths, are necessary to confidently 
determine whether there are temperature variations within the filament itself, and for an 
accurate estimate of the spectral index. 



Kirk et al.l (120071 ) also employed SED fitting to Spitzer and ISO observations to estimate 
the temperatures of numerous cores. We analyze their fluxes in a similar fashion. Though a 
few of the cores did not show definitive temperature drop as shorter wavelength fluxes were 



2 As we discuss in fJSJ a good fit can be found with only three data points if one (or more) of the fit 
parameters (/3, T, or Nkq) is held fixed. 
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excluded, the well studied core B68 shows a clear drop in temperature and an increase in 
the spectral index after short wavelength fluxes are excluded. The fitting results for B68 
is shown in Figure [7(c) |^| The best fit temperatures decreases by ~ 5 K when the 70 and 
90 fim fluxes are excluded in the fit. There is an additional 1 K temperature drop when 
the 160 (Urn flux is excluded. There is also a corresponding increase in (3 from ~ 1.2 to 
~2.4. One interp retation of such a high spectral index is that the dust grains are covered by 
icy mantles (e.g. Kuan et al.lll996l ); this interpretation would be reasonabl e for B68, where 
there is significant molecular depletion onto dust grains (IBergin et al.ll2002l ). The variations 
in the spectral index and temperature estimates when short wavelength fluxes are excluded 
in the fit, along with the relatively small error bars, strongly suggests that there are dust 
temperature (and corresponding inverse density) variations within B68. We thus assign an 
upper limit of 10.8 ± 0.1 K for the coldest region within B68. Our fits also suggest that the 
spectral index of dust in B68 £ 2.4, if that property is constant throughout the core. 



We perform the same test of the data presented by lDupac et al.l (120011 ). They estimated 
the temperature and spectral indices of regions in the Orion complex. One o f the regions 



appea rs to be a dense core without a central source, referred to as "Cloud 2" by lDupac et al. 



(120011 ) . Figure [7J(d) shows our fitting results. For this core, longer wavelength data at 1.2 
and 2.1 mm are available. In this case, the exclusion of short wavelength data reduces the 
temperature by ~ 2 K. The spectral index also increases, from ~2.2 to ~2.5. These results 
are also suggestive of temperature variations within the core, but due to the relatively large 
uncertainties, an isothermal description cannot be ruled out. 

We have found that the systematic exclusion of short wavelength data for dense cores 
results in lower best fit temperatures, and higher best fit spectral indices. A trend of de- 
creasing (3 with increasing T h as been put forward as a physical property of dust grains in 
the ISM by iDupac et al.l (120031 ). They find such a trend by fitting equation ([T]) to observed 
fluxes primarily from PRONAOS/SPM, correspondi ng to wavelengths ( in the range 100 - 
600 yum) near the peak of the SEDs emitted by dust. IDupac et al.l (120031 ) argue that line of 
sight temperature variations only result in a slight variation of (3 with T, and that unreal- 
istically high density contrasts ( ~ 100 x) are neces sary to reproduce the magnitude of the 
inverse T - j3 correlation (see also lDupac et al.ll2002l ). We find that temperature variations in 
realistic cores with a uniform spectral index would show an inverse T - f3 relationship when 
observed at wavelengths near the peak of the emergent SED. Besides resulting in erroneous 
f3 estimates, the best fit parameters vary as short wavelength fluxes are excluded in the fit, 



3 Using all 7 of their fluxes, iKirk et al.l (|2007Q obtain a best fit T=12.5 K; the difference between their 
value and ours arises because they kept j3 at a fixed value of 2. Our results agree when [3 is fixed at that 
value in our fits. 
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with the best fit (3 approaching the correct value when only long wavelengths are considered. 
Further, we also find that modest errors, as low as 5%, in observed fluxes from isothermal 
sources lead to erroneous (3 and T estimates. The trend in the points in a T - (3 plane 
from a number of fits to noisy fluxes also show such an inverse correlation (Fig. [6]). It 
w ould thus be i nform ative to compare the form of the T - (3 relationship we find with that 
of iDupac et al.l (120031 ). We note that though we only consider objects with T < 20 K in this 
work, we investigate the derived inverse T - (3 relationship due to noise in ob servations of 



isoth ermal sources with temperatures up to 100 K in an accompanying paper (jShetty et al. 
2009h . 



3.4-1- A True Inverse T - (3 Correlation? 

To further investigate the derived inverse T - (3 correlations, we again consider the two- 
component medium "2COMPd," with temperatures T\ = 10 K, T 2 = 20 K, and a column 
density ratio N 2 /Ni = 0.1. As shown in FigureEl a modified blackbody SED fit results in a 
T estimate of 23.3 K and a f3 estimate of 0.23, using fluxes measured between A G 10 - 3000 
/im. To compare with realistic observations, we perform the fit to data sets each containing 
fluxes at 5 wavelengths. Two sets of wavelengths are considered: one with fluxes near the 
peak of the SED at 100, 200, 260, 360, and 580 jum, and the other with fluxes in the R-J 
tail of the SED at 850, 1100, 1200, 1500, and 2100 fim. In order to account for the effect of 
a ~5% uncertainty in the observations, due to noise or calibration errors, for example, each 
flux is multiplied by a random value drawn from a Gaussian distribution with a mean of 1.0 
and dispersion of 0.05. We generate 100 sets of fluxes in this manner, and fit equation [1] to 
each set of fluxes, as in the numerical experiments of isothermal sources presented in §3.21 



Figure [8] shows the fitting results from these simulated observations. The fits to the 
observed fluxes at long wavelengths show a scatter in the T - [3 relation, due to the uncer- 
tainties in the fluxes. The fits at shorter wavelengths do not show as much scatter, but the 
best fit 13 and T are very poor estimates of the true (3 or the column temperature of 10.9 K. 
As previously discussed, this inaccuracy at short wavelengths results because the emergent 
SED from a non-isothermal source is not well fit by a single power law modified blackbody 
spectrum. 



Pupae et al.l (120031 ) observe a variety of sources, including dense cores as well as much 
warmer sources. After fitting modified blackbody SEDs to the observed fluxes, with wave- 
lengths < 600 /im, they find that a hyperbolic form in f3(T) represents the fits well. They 
investigated the effect of noise in a model with a range of uncorrelated T - f3 pairs, and 
concluded that such a model is inconsistent with the observed data. The shape of the T - f3 
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correlation in Figu re M from a model with a single (3 is remarkably similar to that shown in 
Pupae et all (120031 ). This suggests that an inverse, and possibly even hyperbolic-shaped, T 
- (3 relationship is not necessarily due to real variations in the dust spectral index with dust 
temperature. The relationship may simply be due to temperature variations along the line 
of sight, along w i th un certainties in the observed fluxes. At the short wavelengths considered 
by iDupac et al.1 (120031 ) , for the warmer sources the observed wavelengths may indeed fall in 
the R-J part of the spectrum. For these sources, the "short wavelength" fits will be analogous 
to the "long wavelength" set in Figure [HI For sources that are dense cores, if they are not 
isothermal, the SEDs at those short wavelengths are not well fit by equation (TjQ); a fit would 
produce a large T estimate, relative to the column temperature, and would underestimate 
(3. 

For isothermal sources ( §3.21) . the addition of noise to the observed fluxes further de- 
grades the parameter estimates. Our analysis indicates that the SED fits to fluxes with A < 
600 /im from various warm isoth ermal sources with T £ 60 K may all produce similar T and 
(3 estimates Jshettv et alJl2009h . Uncertainties in the observed fluxes of starless cores ma y 
be responsible for some of the scatter in the T - (3 diagram shown bv lDupac et al.l ( 2003). 



Howev er, all of the spread is likely not a consequence solely of noise, since IDupac et al. 



( 120031 ) observe a variety of sources. One possibility is that the spectral index varies within 
a source, which is a situation we do not model (see Tabled]). The emergent SED from such 
sources will of course be more complicated, for which alternative analysis techniques may 
provide better parameter estimates. We note that the points in Figure [S] would simply be 
systematicall y offset had we con sidered a source with a different (constant) spectral index 
from [3 = 2 (jShetty et al.l 120091 ) . Thus, had we included multiple sources with different 
(constant) values of (3, the T - (3 diagram would be further populated. Though we cannot 
exclude the possibility that the spectral index of dust decreases with increasing temperature, 
we have shown that simple power-law-modified-blackbody fits to observed data can result in 
misleading T - (3 relationships which appear like those sometimes claimed to be of physical 
origin. 

We have demonstrated that the assumption of isothermality can lead to significant errors 
in estimates of T and /3 from SED fits. Noise contributes to additional uncertainty in the 
estimated parameters. In the next section, we describe the effect of noise on a different 
method commonly used to estimate T and (3, by means of flux ratios, focusing on isothermal 
sources. We then compare the flux ratio estimates to those derived from SED fitting for 
sources with line of sight temperature variations. 
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Flux Ratios 



An alternative to full SED fitting for estimating the temperature is through the use of 
ratios of observed fluxes. For a given observation, the known quantities in equation ([1]) are 
the flux density S u and the beam size Q. With two observations (smoothed to a uniform 
resolution) at different frequencies ^iandz/ 2 , corresponding to wavelengths AiandA 2 , the 
ratio of S Vl /S V2 produces an equation where the two unknown quantities are T and (3: 



"2 



A2 



3+/3 



exp(A T /A 2 ) 
exp(A T /Ai) 



(3) 



where At = hc/kT. The main assumptions made to derive this equation are that T and j3 are 
constant along the line of sight. If (3 is known a priori (or otherwise assumed to be known), 
then equation (j3)) c a n be used to estimate the t emperature from only two ob s ervations (e.g. 



Kramer et all 120031 : iSchnee fc Goodman! 120051 : IWard-Thompson et al.l 120021 : iSchlegel et al. 
1998h . 



Including a third observation at frequency z/3, corresponding to wavelength A3, taking 
ratios using all three flux densities produces 



log 



exp(A T /A 2 ) - 1 
exp(A r /Ai) - 1 



iog (^ 



log 



exp(A T /A 3 



1 



exp(Ar/A 2 



(4a) 
(4b) 



The advantage of using this equation to estimate the temperature is that no assumption for 
the value of (3 is required, though (3 is assumed to be constant along the line of sight. A 
similar equation can be derived for observations at four wavelengths. However, as we discuss 
in £j6l in that case a direct fit of equation ([1]) is reliable. We will hereafter refer to the 
left hand side of equation (j3J) as the "flux ratio," and the right hand side as the "analytic 
prediction." 



Schnee et al.l (120071 ) used fluxes from observations of the starless core TMC-1C at 450, 
850, and 1200 /im in equation (HI) . They found that the errors in the observations would 
have to be 2% in order to accurately estimate the temperature (of an isothermal source). 
The goal of that study, besides mapping the temperature, was also to map the spectral index 
and column density. Once an estimate for the temperature is obtained, the spectral index 
can be estimated through the use of the ratio of any two fluxes as: 



(3 = log 



S V1 exp(A T /Ai) - 1 
S V2 exp(A T /A 2 ) - 1 



/log(^ 



(5) 



-18- 



In principle one can also estim ate the column den sity N, modulo kq, once T and (3 are 
derived from equations (JH)-© (e.g. ISchnee et al.ll2006l ). As discussed in ^J3j however, addi- 



tional assumptions for the opacity and/or the dust-to-gas ratio may be required if extinction 
observations are unavailable. Since our focus is on the temperature and spectral index, we 
do not consider those assumptions here. 



4.1. Determining Temperatures from Two Fluxes 

We begin our analysis of the flux ratio method by considering equation (jHJ) to estimate 
the temperature, given two fluxes at different wavelengths. We consider both an isothermal 
source and a two-component source (as in §3.1.11) . We then compare the fluxes at two 
chosen wavelengths to the analytical prediction from the right-hand-side of equation (J3]). 
We analyze the effect of noise, as well as of an incorrect assumption of j3, on the temperature 
estimate. For the analysis including noise, we add a random component drawn from a 
Gaussian distribution with a chosen dispersion to the flux, and then use these "noisy" fluxes 
in equation Q. We repeat this simple experiment 10,000 times to obtain better statistics 
on the temperature estimates. 

Table [2] shows the temperature estimates using observations with different noise levels, 
and assuming different values of (3. Column (1) shows the wavelengths of the two fluxes. 
Column (2) shows the assumed value of (3 in equation (j3j); the spectral index of the model 
source is 2.0. Column (3) indicates the level of noise added to the fluxes. Column (4) shows 
the true column temperature; for the isothermal source, the column temperature is just the 
actual source temperature. The last column gives the derived temperature. For the fluxes 
that are altered by noise, we show the la distribution in the estimated temperatures. 

For an isothermal source, the shorter wavelength pair is less sensitive to noise and/or 
errors in the assumed value of f3. For the two-component source "2COMPd" (see Fig. 
[3]) without noise, the longer wavelength pair gives a (slightly) more accurate measure of 
the column temperature than the shorter wavelength pair, when the correct value of (3 
is assumed. However, when modest levels of noise are included in the fluxes, the long 
wavelength pair produces temperature estimates that deviates from the true temperature 
(or column temperature) more than the short wavelength pair. Additionally, the longer 
wavelength pair is more sensitive to the assumed value of (3. 

These trends are as expected given the shape of the modified blackbody (Fig. [I]). As 
we demonstrated in fusing a direct SED fit, fluxes in the R-J part of the spectrum provide 
more accurate T and (3 estimates, though the fits are more sensitive to noise. Similarly, 



-19- 



the flux ratio involving wavelengths in the R-J part of the spectrum is more sensitive to 
the assumed value of the spectral index, as well as noise, indicating that short wavelength 
fluxes are preferable. However, when large temperature gradients are present (see §3.1.1l and 
§3.1.2p . the temperature estimate from the short wavelength pair will deviate significantly 
from the column temperature. 

For isothermal sources, the source temperature, as well as the assumed value of (3 
determines how well fluxes at different wavelengths could recover the temperature. For noise 
free observations of a 10 K source at 450 and 850 /im, and when {3 is assumed to be within 
5% of the source value, T can be recovered within 10%. The accuracy of the T estimate 
degrades to ~20% with 1200 and 2100 /im observations. For warmer sources (T >> 10 K), 
450 and 850 /im fluxes would be well into the R-J part of the spectrum. These fluxes would 
thus be more sensitive to noise and (3 than those 450 and 850 /im fluxes from a 10 K source. 
Lower temperature sources always give better (i.e. less uncertain) temperature estimates. 
The flux ratio method, however, still requires a reasonable assumption for (3. 

For non-isothermal sources, the estimates become more uncertain. As shown in Table 
HI low noise levels and an accurate assumption of (3 can reasonably recover the column 
temperature of a simple 2 component medium. Since equation ([3]) is derived from the 
isothermal assumption, the estimated temperature becomes less accurate for more complex 
sources, even when comparing with the column temperature. In our subsequent analysis of 
flux ratios, we will hereafter concentrate on deriving temperatures of isothermal sources. 



4.2. Determining Temperatures from Three Fluxes 

We next investigate the accuracy in estimating temperatures from observations at 3 
wavelengths, using Equation (j3J). We first consid er observatio n s at 450, 850, and 1200 



/im, the three "popular" wavelengths used in the ISchnee et all (120071 ) study. Figure (jHj) 
shows the analytic prediction used to determine T from expression (14b j) for temperatures 
between 1 - 100 K. At temperatures T £ 7 - 10 K, the analytic prediction is very sensitive 
to the temperature. Thus, even though errors in the fluxes will produce an inaccurate 
value in expression (I4"a|) for comparison with the analytic prediction of expression (14b p . the 
derived temperature will still be close to the actual temperature of the emitting medium. 
At temperatures T £ 7 - 10 K, however, the analytic prediction is not very sensitive to the 
temperature. Thus, even small errors in the flux, due to noise and/or other observational 
uncertainties, will result in grossly erroneous temperature estimates. 



To investigate the effect of noise on the determination of temperature using equation 
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IU we have run a number of Monte-Carlo simulations. In these simulations, the emergent 
flux of a source at constant temperature (eqn. [1]) is modified by some chosen level of 
Gaussian noise, representing random errors in real observations. These "observed" fluxes 
at three wavelengths are used in expression (|4ap to compare with the analytic prediction of 
expression (jib]) . 



The uncertainties in the observations of lSchnee et al.l (120071 ) were estimated at 12%, 4%, 



and 10% for the 450, 850, and 1200 /im observations, respectively. We use those uncertainties 
in our first simulated observations. We "observe" the extended source at 10,000 positions 
(or, equivalently, 10,000 times at a single location on the sky) at those three wavelengths. 
The flux ratios, using the noisy fluxes in expression (j4aj) . from two sources at 10 K and 5 
K are shown in Figure (TDJ The mean values of the flux ratio recovers the true temperature 
reasonably accurately. Also marked are the ±2<r levels (as well as the ±3<r levels for the 
5 K source). At 10 K, even at the +2<r level, the flux ratio does not lie in the 1 - 100 
K range of the analytic prediction shown in Figure [9j However, for the source at 5 K, at 
the 3a level the derived temperature only differs from the true temperature by a factor of 
~ 2.4. At 10 K (and higher) temperatures, Figure [9] shows that the analytic prediction 
does not vary much with temperature, so any error in the flux ratio will correspond to a 
temperature that deviates significantly from the true temperature. For a 5 K source, the 
analytic prediction varies significantly with slight variations in temperature, so errors in 
the flux ratio will still produce reasonably accurate temperature estimates. This test has 
shown that with observations at 450, 850, and 1200 /im, one can only be confident in the 
ratio method if the estimated temperatures are ^ 5 K. For temperatures greater than the 
turn-over temperature in Figure [9] (~ 7 - 10 K), the derived temperatures cannot be deemed 
accurate. 

The temperature at which analytic prediction shifts from a strong temperature depen- 
dence to a weak temperature dependence is determined by the particular wavelengths used 
in expression (I4bl) . In this sense, the ideal set of three wavelengths would shift the turn- 
over temperature to much higher values. In order to test the sensitivity of the turn-over 
temperature to the particular values of wavelengths, we have run a series of Monte-Carlo 
simulations as described above; we varied both the set of 3 wavelengths, as well as the 
(constant) temperature of the source. 

In order to locate the turn-over temperature, we define a threshold temperature Tth 
such that a ±3<r range in the flux ratios corresponds to estimated temperatures T est that are 
within a factor of 2 of the actual source temperature T . For example, consider a source at T 
= 21 K observed at three given wavelengths. If the range of derived temperatures included 
in the ±3cx level of the flux ratio includes temperatures > 42 K, then we know that T t h < 21 
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K for that set of wavelengths. If we then run another simulation on a source with To = 20 
K, with the same three wavelengths, and find that the maximum derived temperature in 
the 3a range is T est < 40 K, then we set T t ^ = 20 K for this set of three wavelengths. Our 
definition of T t h is arbitrary, and can of course be set to correspond to a higher accuracy 
temperature estimate; the goal here is simply to compare how well observations with varying 
wavelengths can reproduce source temperatures to a chosen level of accuracy. We note that 
at the -3cr level, the difference of \T est — T \ is always less than that difference at the +3cr 
level, due to the logarithmic functional form of the analytic prediction; at lower flux ratios, 
corresponding to lower values of the analytic prediction (eqn. jib]), the derived temperature 
is less sensitive to uncertainties in the fluxes (see Figs. 191-fTOl). 

Figure [TT] shows the threshold temperatures, for given wavelength ratios A3/A2 and 
A 2 /Ai, where Ai < A 2 < A 3 @ For simplicity, we assume a 10% noise level in the fluxes at all 
three wavelengths. A clear trend is immediately apparent in Figure [TTJ The highest values 
of T t h occur when the two ratios A3/A2 and A2/A1 are simultaneously large. However, for a 
given ratio A3/A2 or A2/A1 there is a limit to how large the other ratio can be, beyond which 
T t h decreases. From Table El the highest T t h is achieved when Ai = 70 jum and A 3 > 2000 /im. 
At those wavelengths, both the Wien and R-J regimes of the SED for sources with T £ 75 
K are sampled. 

Given two wavelengths that sample the Wien and R-J limits, the maximum temperature 
that can be reliably found is set by the middle wavelength A2: the maximum temperature is 
roughly that determined by Wien's displacement law using A2 as the wavelength correspond- 
ing to the flux peak. This general trend evidently dissolves as one or more of the boundary 
wavelengths Ai or A3 approaches A2. Since the wavelength corresponding to the SED peak is 
inversely proportional to the temperature, an increase of all three wavelengths by a constant 
factor would result in a decrease in T t h by that same factor. Table El indeed shows such a 
trend. In short, the exact value of T t h is dependent on all three wavelengths, as would be 
expected. 

We have shown that given three wavelengths, one could determine a threshold tem- 
perature above which the ratio method will not be able to accurately derive the source 
temperature (to some chosen level of accuracy), due to uncertainties in the observations. In 
general, lower temperatures regions, such as cold, dense cores, can be more accurately mea- 
sured through this method. In principle one could also determine which three wavelengths, 



4 Tablc[3]in the Appendix explicitly lists the threshold temperatures T t h for various sets of wavelengths, 
along with the wavelength ratios. Our choice of the three wavelengths span the range 70 - 3000 fim. These 
particular wavelengths are shown because many of them arc included in the wavebands of the upcoming 
Herschel and Planck missions. 
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and the associated uncertainties, are optimal to estimate a given temperature. However, 
such an analysis would not be as practical, since observers do not have an arbitrary choice 
as to what wavelengths they can observe, nor to any desired level of accuracy. After a brief 
discussion on estimating ft, we will show that fitting the SED directly to estimate the tem- 
perature is much more accurate than the ratio method, regardless of what three wavelengths 
are observed. 



4.3. Estimating ft 

In the ratio method, once a temperature estimate is obtained, it may be used in equation 
(J5J) to estimate ft. Only two fluxes are required in equation ((Sj). For isothermal sources, any 
two fluxes will give good estimates of ft. Since the SED is more sensitive to noise at longer 
wavelengths, shorter wavelength fluxes will produce more accurate ft estimates. However, 
determining whether a source is isothermal itself is not trivial; as discussed in §3.41 with 
numerous fluxes this can be accomplished by systematically excluding short wavelength 
fluxes in a direct SED fit. When a source is not isothermal, the resulting temperature 
estimates from short wavelength fluxes will be highly inaccurate. The use of these incorrect 
temperatures in equation (jSJ) will also produce incorrect ft estimates. Thus, an accurate 
estimate of T is required before accurately estimating ft through equation (j5J). 

Uncertainties in observations in the R-J tail will result in highly inaccurate ft estimates 
using the flux ratio method. For example, consider an isothermal source that is observed at 
three wavelengths, 450, 850, and 1200 //m. Fluxes that are inaccurate by a mere 3% can 
produce ft estimates that deviate from the actual value by 25%. For the simple 2 component 
source considered in §3.1.11 (and §4.10 . and with flux uncertainties of only 5%, the estimated 
ft can be inaccurate by ~50%. Thus, the flux ratio method gives highly uncertain estimates 
of ft. 

5. Comparison of Fitting and Flux Ratio Methods involving Three Fluxes 

In this section, using simulated observations of a dense core with temperature and 
density gradients, we compare the temperature estimate from direct SED fit to that derived 
from a flux ratio. Fluxes are only "observed" at three wavelengths, and include uncertainty 
due to noise, requiring that one of the three parameters in the fit (T, ft, or Nkq) is held 
fixed. 

One approach to carry out the comparison would be to construct SEDs throughout the 
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volume of the core, using equations ([I])-©, assuming a fixed value for (3. We could then 
integrate along all lines of sight to obtain the emergent intensity at three wavelengths, and 
then carry out the fitting, as we did in §31 or use the flux ratio method to estimate the 
temperature. 

An alternative approach, which we choose to use here, is to utilize a radiative transfer 
(RT) code. With an RT simulation, the dust properties can be s et using real dust op - 
tical constants. We thus use the radiative transfer code MOCASSIN (jErcolano et al.l 120051 ) . 
MOCASSIN is a 3D code that uses a Monte Carlo ap proach to the transfer of radiation through 
an arbitrary medium (jErcolano et al.ll2003l . 120081 ). 



Our spherically symmetric model core emits radiation according to the emissivity prop- 
erties of the dust. The physical properties of the core, such as the density and temperature 
profiles, as well as the dust properties, are held constant. Since we are only interested in 
the form the observed SED, and not the detailed distribution, the computational grid only 
has 16 x 16 x 16 zones, which follows the propagation of radiation through one-quarter of 
a spherical model core, with appropriate boundary conditions. The overall dimensions of 
the grid is ~0.1 pc 3 ; the resolutio n of an individual zone is c omp arable to the r e soluti on of 
the TMC-1C maps presented by ISchnee &: Goodman! (120051 ) and ISchnee et al.l (120071 ). As 
radiation (or energy packets) from all spatial locations traverses through the cloud, it is 
absorbed and re-emitted by the dust. However, we maintain the original temperature of 
the dust throughout the simulation. In this sense, the temperature of the core is set by 
some external source, such as an ambient interstellar radiation field, which is not explicitly 
included in our model. Radiation that emerges out of the core contributes to the "observed" 
flux. The emergent flux will thus be proportional to the density along the line of sight; the 
resulting observed 2D map of the core will scale with the 3D density integrated over the line 
of sight. We can then apply the flux ratio and fitting methods to estimate the temperature 
at each location on the 2D map. 

We construct a core with a Bonnor-Ebert like density profile, where the density is con- 
stant in the central regions, but then drops off as the square of the radius. The temperature 
of this core is also constant in the central regions, but then increases logarithmically with 
radius. The temperature of the core varies from ~7.5 K at the inner regions, to ^ 12 K at 
the edge of the core. Such values span temperatures near the threshold temperature of the 
flux ratio method using 450, 850, and 1200 fim fluxes, as well as temperatures well into the 
regime where the method becomes extremely sensitive to noise. 



Figure [TJ] shows the core temperatures, both the actual temperature and the recovered 
ones. The ratio method produces a spread of temperatures at any given radius, due to noise 
introduced by the finite number of photons tracked in the simulation. A similar simulation of 
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a core with constant density and constant temperature produces normally distributed fluxes 
with a ~ 10%. Despite the scatter in the observed fluxes, a general trend of decreasing 
temperature with decreasing radius is apparent. Towards the innermost regions, the flux 
ratio method overestimates the temperature more so than in the outer regions. This occurs 
because in the inner regions of the 2D observed map there is still a contribution to the flux 
from the warmer dust at larger radii, from matter lying above the colder central regions. At 
large radii, line of sight variations in temperature and density is minimized, so the mean of 
the flux ratio estimated temperatures corresponds well with the actual temperatures. Using 
the estimated temperatures from the flux ratio method in equation [5] gives (3 = 2.0 ± 0.1. 
This is in good agreement with our choice of silicate grains in the simulation, which should 



have (3=2 flDraine fc Led 119841 ) . 



In performing a fit to only three fluxes at each location, one of the free parameters must 
be held fixed. Since the flux ratio method produces a (3 estimate that is well constrained, even 
though T shows relatively higher levels of scatter, we hold (3 fixed at that value. As shown 
in Figure [121 the temperatures obtained through such a fit are also overestimated at small 
radii. But, the scatter in the estimated temperatures is much smaller using a constrained 
SED fit (i.e. (3 fixed) compared with the flux ratio method. 

Also plotted on Figure [12] is the true column temperature (i.e. density weighted tem- 
perature). Each location (or grid zone) in the 2D map corresponds to a line of sight through 
the 3D core. We integrate the temperature along the line of sight at each location in the 2D 
map, weighted by the 3D densities used to construct the core. Evidently, the fit temperatures 
coincide remarkably well with the column temperature. The range in temperatures from the 
flux ratio method is also spread more evenly about the column temperature than the actual 
3D temperature. At the innermost regions, there remains a slight offset. In general, though, 
the temperature estimates are certainly more representative of the "column temperature" 
than the true temperature of the core. This is not too surprising, because the observed flux 
indeed encodes information about all matter along a line of sight. 

Had we assumed a different (and thus incorrect) value of (3 that was held fixed in the fit, 
the best fit temperature would be systematically offset from the column temperature. For a 
10 K isothermal source, fixing f3 at 1.7 and 2.3 would produce fit temperatures of ~ 12 and 
~ 9 K, respectively, using fluxes at 450, 850, and 1200 /jm. The results of our simple tests 
suggest that though the flux ratio method can give highly uncertain temperature estimates 
from three fluxes, due to line of sight temperature variations, those te mperatures can stil l 



provide a decent first estimate for the (mean) spectral index (as found by lSchnee et al.l (120071 ) 
for isothermal sources) . This estimate for (3 can then be held fixed in a fit which will recover 
temperatures with less scatter, and in close agreement with the column temperature. 
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6. Discussion 

A flux ratio method to estimate the temperatures can also be constructed with four 
wavelengths. The form of the analytical prediction for four wavelengths is similar to ex- 



pression (14b j) . We have also performed Monte Carlo simulations to test whether a flux ratio 
method using four wavelengths produces more accurate estimates of the source temperature, 
compared with a flux ratio method involving only three wavelengths. The general trends 
we found for three wavelengths remains: lower wavelength observations produce less scatter 
in the estimated temperatures, and that the fluxes should sample different regions of the 
spectrum to obtain decent temperature estimates. However, the maximum temperature at 
which the turn over occurs for any set of 4 wavelengths listed in Table [3] is ~ 65 K, for our 
definition of T t h in §4.21 

When fluxes at four wavelengths are available, however, a direct fit may be employed to 
estimate (3 (and the absolute scaling Nk ) along with the temperature. Since determining j3 
through the the flux ratio (eqn. [S]) could give contradictory estimates depending on which 
the fluxes are used, a direct SED fit is preferable. We did not find any advantage of using 
the flux ratio method to determine the temperature using four wavelengths compared with 
a direct fit of a modified blackbody SED. 

In all of our tests, we have only considered sources with constant spectral indices. A line 
of sight may also have a variations in j3, and would further complicate the estimation of dust 
temperature. It may be reasonable to assume that the dust emissivity is constant within a 
core, where temperatures only vary by ~ 10 K. But for lines of sight extending through a 
wider range in density and temperature, such as the lower density, warmer gas surrounding 
sites of recent star formation, assigning a single value for the spectral index may lead to 
errors in determining the temperature. A thorough investigation of spectral index variations 
over a range of environments would be required to quantify its effect on the emergent SED. 



7. Summary 

We have investigated the effect of noise and line of sight temperature variations on 
two common methods used to estimate the dust temperature and spectral index of cold 
star forming cores using continuum observations. One method is a direct fit to a modified 
blackbody spectrum. The second method involves the use of flux ratios. 

We demonstrate that employing an isothermal modified blackbody equation (eqn. [T]- 
[2]) may lead to highly inaccurate dust temperature and spectral index estimates. Least 
squares SED fits to fluxes in the R-J regime, as opposed to the Wien regime, may provide 
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accurate spectral index and density weighted temperature, or column temperature, estimates. 
For conditions typical of starless cores, fluxes in the R-J regime have wavelengths ^ 600 
/im. However, the fits to fluxes in the R-J regime are rather sensitive to observational 
uncertainties, such as noise. 

The flux ratio method may also provide inaccurate parameter estimates due to line of 
sight temperature variations, and is also very sensitive to noise. In a comparison of the flux 
ratio and least squares fitting methods when only three fluxes are available, we find that a 
direct fit with the spectral index held fixed provides more accurate estimates of the column 
temperature. The flux ratio method can be initially used to estimate the value of the spectral 
index to be held fixed for a least-squares SED fit. 

We summarize our main findings in more detail here: 

1) Line of sight temperature variations can lead to inaccurate temperature and spectral 
index estimates when fitting a power-law-modified blackbody SED to observed fluxes. Near 
the SED peak of sources with temperature variations, the spectrum is poorly fit by an 
isothermal spectrum. For longer wavelength observations in the Rayleigh- Jeans regime of 
the spectrum, and with minimal observational uncertainties, a fit can accurately recover the 
spectral index (if it is constant), and provides a good estimate of the upper limit of the 
column temperature. However, at these long (Rayleigh- Jeans) wavelengths a fit is extremely 
sensitive to noise. 

2) Short wavelength observations (A ^ 600 /im) are still useful, for they can indicate 
whether an observed source contains temperature variations. For starless- core- like sources 
with temperature variations, the resulting fit T decreases and fit 3 increases when system- 
atically exclu ding short wav e lengt h flu xes from the fit. Pub lished data of sources in Taurus 
and Orion by IStepnik et al.l (120031 ) and lDupac et al.l (120011 ). respectively, show these appar- 
ent trends, but an isothermal description with no systematic var iations in 8 still c annot be 
strictly ruled out, due to the uncertainties. Observed fluxes by Kirk et al.l (120071 ) of B68, 
though, produce lower fit temperatures and higher fit spectral indices when short wavelength 
fluxes are omitted, strongly suggestive of dust temperature variations along the line-of-sight. 
We estimate an upper limit of 10.8 ± 0.1 K for the temperature of the coldest region within 
B68; and if the spectral index is constant throughout the core, then we estimate 3 > 2.4. 

3) SED fits to fluxes in the Rayleigh- Jeans regime are very sensitive to noise, even for 
isothermal sources. The fits may produ ce a spurious inverse T-3 relationship, similar to the 
trend discussed by iDupac et al.l (120031 ). SED fits may be more accurate when fluxes with 
wavelengths that span the SED peak are available, compared with fits to fluxes solely in 
the Rayleigh- Jeans regime. However, fits to fluxes near the SED peak would be inaccurate 



-27- 



if the source contains line of sight temperature variations. In general, for objects that are 
cool pockets in higher density regions, such as starless cores, SED fits that produce higher 
temperatures also (artificially) give lower spectral indices. 

4) We find that, due to noise uncertainties in any observation, the flux ratio method 
is most accurate for emission originating from cold isothermal regions. For a source with a 
constant temperature, there may be a range in the estimated temperature, due to the uncer- 
tainties in the observations. At low temperatures, this spread is small; at higher temperatures 
the spread can be rather large, rendering the temperature estimate highly inaccurate. The 
precise temperature, or threshold temperature, for which the method shifts from relatively 
accurate to inaccurate is dependent on the observed wavelengths (as well as the desired level 
of accuracy, see Fig. [LTJand Table ED- For example, for fluxes at 450, 850, and 1200 /im, the 
flux ratio method can provide accurate temperature estimates only for sources with T £ 7 
- 10 K (Fig. ED. 

5) Using Monte Carlo simulations, we quantified the dependence of the turn over tem- 
perature on the set of observed wavelengths. Ideally, as one might intuitively expect, two of 
the three wavelengths should sample the Wien and the Rayleigh- Jeans regime of the SED, 
with the final wavelength lying at intermediate values. Further, a greater separation be- 
tween the wavelengths results in more accurate temperature estimates. In general, higher 
temperatures can be more accurately measured when the observations include short wave- 
length far-infrared observations A ^100 /im; at short wavelengths, however, there may be 
a contribution from stars and transiently heated very small grains to the observed flux. 

6) A reasonably accurate estimate of j3 can be obtained from the mean of the estimates 
derived from the flux ratio method involving three fluxes. For fluxes at 450, 850, and 1200 
/im, with ~10% uncertainties, (3 can be estimated to within 5% of the true source value. 
This value can then be held fixed in a constrained SED fit to the three fluxes to estimate the 
temperature with less scatter than an estimate from the flux ratio method (Fig. [T2|) . With 
four or more observations j3 may be one of the free parameters in the fit (and, of course, the 
fit is better constrained). 

7) The temperatures estimated through the SED fit and ratio methods, however, cannot 
be used to assign the absolute temperature to a given 3D location in a cold core. The 
projected SED contains information from all emitting matter along any line of sight. The 
measured temperature is more representative of the column temperature. In this regard, the 
estimated temperatures provide an upper limit for the coldest temperature along the line of 
sight. 
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A. Appendix 

Table [3] shows the threshold temperatures T t h given fluxes at three wavelengths (see § 
14. 2p . The ratios A3/A2 and A2/A1 are also provided. The table explicitly shows the quantities 
used to produce Figure [HQ 
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Table I. Overview of §3 - §6 



Scenario 


§3.1 


§3.2 


§3.3 


§3.4 


§4.1 


§4.2 


§4.3 


§5 


§6 


Least Squares Fitting 
Flux Ratio Method 


/ 


/ 


/ 


App a 


/ 


/ 


/ 


/ 
/ 




T constant 
Two T Medium 
Gradient in T 


/ 
/ 


/ 


/ 




/ 
/ 


/ 




/ 




(3 constant 
(3 variable 


/ 


/ 


/ 


App a Fixed b N/A c Derived d 


/ 


/ 


Two Fluxes 
Three Fluxes 

# Fluxes > 3 

# Fluxes > 3 


/ 


/ 


/ 


App a 


/ 


/ 




/ 


/ 


Without Noise 
With Noise 


/ 


/ 


/ 


App a 


/ 
/ 


/ 




/ 





Application: SED fits to fluxes from published observations 

b Choice of (3 required in method 

c f3 is not required or recovered from method 

d /3 can derived from flux ratio method involving 2 fluxes 
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Table 2. Temperature Estimates from 2-Flux Ratio Method 



Observed Ai, A 2 


Assumed f3 & 


O'obs 


Source T cot 


Derived T 


(H 




(%) 


(K) 


(K) 


Isothermal: 










450,850 


1.9 





10 


10.5 


450,850 


2.0 





10 


10.0 


450,850 


2.1 





10 


9.5 


450,850 


1.9 


10 


10 


10.6 ± 1.3 


450,850 


2.0 


10 


10 


10.0 ± 1.2 


450,850 


2.1 


10 


10 


9.6 ± 1.0 


1200,2100 


1.9 





10 


12.1 


1200,2100 


2.0 





10 


10.0 


1200,2100 


2.1 





10 


8.6 


1 9nn 9i nn 


1 Q 


1 n 


1 n 


14 R 4- 1 9 fi 


1200,2100 


2.0 


10 


10 


12.1 ± 9.4 


1200,2100 


2.1 


10 


10 


10.2 ± 7.0 


2-Component Source (2COMPd): 










450,850 


1.9 





10.9 


12.5 


450,850 


2.0 





10.9 


11.7 


450,850 


2.1 





10.9 


11.0 


450,850 


1.9 


10 


10.9 


12.7 ± 2.0 


450,850 


2.0 


10 


10.9 


11.8 ± 1.7 


450,850 


2.1 


10 


10.9 


11.1 ± 1.5 


1200,2100 


1.9 





10.9 


13.8 


1200,2100 


2.0 





10.9 


11.1 


1200,2100 


2.1 





10.9 


9.3 


1200,2100 


1.9 


10 


10.9 


15.4 ± 13.4 


1200,2100 


2.0 


10 


10.9 


13.5 ± 10.9 


1200,2100 


2.1 


10 


10.9 


11.3 ± 8.2 



a Model Spectral index (3 



Table 3. Threshold Temperatures in Ratio Method 





\ 2 


Aq 


A3/A2 


Ao/Ai 


± th 


70 


110 


170 


1.55 


1.57 


18 


70 


110 


350 


3.18 


1.57 


38 


70 


110 


450 


4.09 


1.57 


41 


70 


110 


850 


7.73 


1.57 


47 


70 


110 


1200 


10.91 


1.57 


49 


70 


110 


1380 


12.55 


1.57 


50 


70 


110 


2100 


19.09 


1.57 


51 


70 


110 


3000 


27.27 


1.57 


52 


70 


170 


350 


2.06 


2.43 


39 


70 


170 


450 


2.65 


2.43 


47 


70 


170 


850 


5.00 


2.43 


55 


70 


170 


1200 


7.06 


2.43 


57 


70 


170 


1380 


8.12 


2.43 


58 


70 


170 


2100 


12.35 


2.43 


59 


70 


170 


3000 


17.65 


2.43 


60 


70 


350 


450 


1.29 


5.00 


18 


70 


350 


850 


2.43 


5.00 


47 


70 


350 


1200 


3.43 


5.00 


54 


70 


350 


1380 


3.94 


5.00 


55 


70 


350 


2100 


6.00 


5.00 


59 


70 


350 


3000 


8.57 


5.00 


61 


70 


450 


850 


1.89 


6.43 


37 


70 


450 


1200 


2.67 


6.43 


50 


70 


450 


1380 


3.07 


6.43 


52 


70 


450 


2100 


4.67 


6.43 


57 


70 


450 


3000 


6.67 


6.43 


60 


70 


850 


1200 


1.41 


12.14 


22 


70 


850 


1380 


1.62 


12.14 


28 


70 


850 


2100 


2.47 


12.14 


44 


70 


850 


3000 


3.53 


12.14 


52 


70 


1200 


1380 


1.15 


17.14 


9 



Table 3 — Continued 



Ai 


A 2 


A 3 


/ Ao 


A 2 /Ai 


T t h 

Lib 


70 


1200 


2100 


1.75 


17.14 


29 


70 


1200 


3000 


2.50 


17.14 


42 


70 


1380 


2100 


1.52 


19.71 


23 


70 


1380 


3000 


2.17 


19.71 


37 


70 


2100 


3000 


1.43 


30.00 


19 


110 


170 


350 


2.06 


1.55 


18 


110 


170 


450 


2.65 


1.55 


21 


110 


170 


850 


5.00 


1.55 


27 


110 


170 


1200 


7.06 


1.55 


29 


110 


170 


1380 


8.12 


1.55 


30 


110 


170 


2100 


12.35 


1.55 


31 


110 


170 


3000 


17.65 


1.55 


33 


110 


350 


450 


1.29 


3.18 


11 


110 


350 


850 


2.43 


3.18 


30 


110 


350 


1200 


3.43 


3.18 


36 


110 


350 


1380 


3.94 


3.18 


37 


110 


350 


2100 


6.00 


3.18 


42 


110 


350 


3000 


8.57 


3.18 


45 


110 


450 


850 


1.89 


4.09 


24 


110 


450 


1200 


2.67 


4.09 


32 


110 


450 


1380 


3.07 


4.09 


34 


110 


450 


2100 


4.67 


4.09 


41 


110 


450 


3000 


6.67 


4.09 


44 


110 


850 


1200 


1.41 


7.73 


14 


110 


850 


1380 


1.62 


7.73 


19 


110 


850 


2100 


2.47 


7.73 


29 


110 


850 


3000 


3.53 


7.73 


36 


110 


1200 


1380 


1.15 


10.91 


6 


110 


1200 


2100 


1.75 


10.91 


20 


110 


1200 


3000 


2.50 


10.91 


29 


110 


1380 


2100 


1.52 


12.55 


15 
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Table 3 — Continued 



Ai 


A 2 


A 3 


A3/ A2 


A 2 /Ax 


T a 
J-th 


110 


1380 


3000 


2.17 


12.55 


25 


110 


2100 


3000 


1.43 


19.09 


13 


170 


350 


450 


1.29 


2.06 


6 


170 


350 


850 


2.43 


2.06 


17 


170 


350 


1200 


3.43 


2.06 


20 


170 


350 


1380 


3.94 


2.06 


21 


170 


350 


2100 


6.00 


2.06 


23 


170 


350 


3000 


8.57 


2.06 


25 


170 


450 


850 


1.89 


2.65 


15 


170 


450 


1200 


2.67 


2.65 


19 


170 


450 


1380 


3.07 


2.65 


21 


170 


450 


2100 


4.67 


2.65 


24 


170 


450 


3000 


6.67 


2.65 


27 


170 


850 


1200 


1.41 


5.00 


10 


170 


850 


1380 


1.62 


5.00 


12 


170 


850 


2100 


2.47 


5.00 


19 


170 


850 


onnn 

3000 


3.53 


5.00 




170 


1200 


2100 


1.75 


7.06 


14 


170 


1200 


3000 


2.50 


7.06 
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a T t /j chosen as the temperature at which ±3<r is 
within a factor of 2 of the source temperature. 
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500 200 600 \WQ 

Fig. 1. — Emissivity modified blackbodies with different spectral indices (3, but constant 
column density, from a 20 K source. Dashed SED is from a 5 K source with (3 = 2. Thin 
solid vertical lines indicate the peak wavelength of the SEDs. 
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Fig. 2. — Modified blackbody SEDs from a two-component source with temperatures T— 10 
K (solid) and T = 15 K (dashed). The column density of the cooler source is a factor of 10 
larger than that of the warm source. The boundary of the shaded region is the integrated 
SED, from a line of sight containing both sources. 
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Fig. 3. — Best fit T (solid) and (3 (dashed) to emergent SEDs from two- component sources, 
using fluxes in different wavelength ranges. The top and middle rows show the best fit T 
and /?; the left ordinate shows T and the right ordinate shows (3. The dark line shows the 
column temperature T co i, and the light line indicates the /3. The bottom row shows the 
corresponding wavelength ranges of the fluxes used in each fit. The emergent SED from 
"2COMPb" is shown in Figure [2J 
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Fig. 4. — Best fit T (solid) and (3 (dashed) for two cores with and without an envelope, 
using fluxes in different wavelength ranges, as in Figure [31 For Core 1, T varies between 
8-12 K and N between 2xl0 21 - 1.25 x 10 22 cm" 2 ; For Core 2, T varies between 5-12 K and 
N between 2xl0 21 - 1 x 10 23 cm" 2 . For the Envelope, T = 20 K and N = 10 21 cm' 2 . The 
dark line shows the column temperature, and the light line indicates the spectral index. 
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Fig. 5. — Actual and fit SEDs from model Core 2 (see text). The boundary of the shaded 
region is the emergent SED of the core. The dashed SED shows the best fit to fluxes between 
100 - 600 fim (marked by squares). The solid line shows the best fit to fluxes between 1000 - 
3000 /iin (marked by triangles). The green and red lines marks the extent of the wavelength 
ranges used in the two fits. 



2.6 



2.4 



2.2 



~i — i — i — i 



A 
A 



2.0 



% 



A 



~i 1 1 1 1 1 r~ 



□ 
□ 



_i i i i_ 



□ 

□ 

□ □ 



mfco Jog 



— i — i — i — i 

a 10 K 
□ 20 K 



-i i | i i i r- 

100-600 fj,m 

500-1000 /jm 
1000-1500 fj.m 



~~\ 1 1 r~ 




A 

A A 



_J I I I I I L_ 



□ D D □ □ □ 



□ □ 



A 

A A 
A 



_l l l l_ 



20 



A I 

30 



_l I I l_ 



Fit T (K) 



40 



Fig. 6. — Best fit T and (3 from Monte Carlo simulations of noisy fluxes from 10 K (tri- 
angles) and 20 K (squares) isothermal sources. The vertical lines indicate the true source 
temperatures, and the horizontal line marks the true spectral index. Different wavelengths 
fluxes were considered in each fit: 100-600 /im (blue), 500-1000 /im (green), and 1000-1500 
/im (red). Gaussian distributed noise is added to each flux, with a = 5%. 
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Fig. 7. — Best fit T (top row) and /3 (middle row) to small number of observed fluxes, marked 
by the line segments in the bottom row. a) A 2 component source (with T\ = 10 K, T2 = 
15 K, and N 2 /N 1 =0.02), observed at A = 60, 100, 200, 260, 360, and 580 ^m. The bottom 
panel also shows the emergent SED from this source, with the abscissa corresponding to 
\ok(S„); b) F il amen t in Taurus, observed at A = 60, 100, 200, 260, 360, and 580 /im by 
Stepnik et all J2003h : c) B68, observed at A = 70, 90, 160, 170, 200, 450, and 850 /im by 
Kirk et all (|2007h: a nd d) Core in Orion, observed at A = 200, 260, 360, 580, 1200, and 2100 



/jm by iDupac et al.l (120011 ) . 
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Fig. 8. — Best fit (3 and T to observed fluxes from the 2 component source "2COMPd" (with 
7i = 10 K, T 2 = 20 K, and N 2 /N 1 =0.1; see Fig. [3]). Five fluxes are used in each fit: 850, 
1100, 1200, 1500, and 2100 fim (squares), or 100, 200, 260, 360, and 580 ^m (triangles). 
Each flux includes a small (Gaussian distributed) random component, with a = 5%. The 
lines indicate the model parameters (3 = 2 and column temperature = 10.9 K. 
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Fig. 10. — Histogram of flux ratios (LHS of eqn. [1]), including noise, for sources at (a) 10 
K and (b) 5 K. Noise levels were set at 12%, 4%, and 10% for the 450, 850, and 1200 /an 
observations, respectively. Mean flux ratio corresponds to (left) 9.8 K, and (right) 5.0 K. 
Lines on the 10 K histogram show the ±2<r level, corresponding to temperatures of > 100 K 
and 5.5 K (see Fig. [9]). Lines on the 5 K histogram show ±2(7 levels, corresponding to 8.0 
K and 3.6 K, and ±3cr, corresponding to 11.9 K and 3.1 K. 
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Fig. 11. — Threshold temperature T t h (see §4.21) from the flux ratio method, for given ratios 
A 3 /A 2 and A 2 /Ai, where Ai < A 2 < A 3 . 
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Fig. 12. — Comparison of flux ratio and SED fitting methods to estimate fluxes. Solid black 
line shows actual (3D) temperature profile of core. The red squares show temperatures 
derived using flux ratio method. Blue crosses show best fit SED temperatures assuming a 
fixed value of (3 = 2.0. The green line shows the integrated temperature along the line of 
sight, weighted by the density, or "column temperature." 



